A matched expansion approach to practical 
self-force calculations 



Warren G. Anderson and Alan G. Wiseman 

Department of Physics, University of Wisconsin - Milwaukee, P.O. Box 413, 
Milwaukee, Wisconsin, 53201, U.S.A. 

E-mail: warrenSgravity . phys . uwm . edu , agwOgr avity . phys . uwm . edu 

Abstract. We discuss a practical method to compute the self- force on a particle 
moving through a curved spacetime. This method involves two expansions to calculate 
the self-force, one arising from the particle's immediate past and the other from the 
more distant past. The expansion in the immediate past is a covariant Taylor series 
and can be carried out for all geometries. The more distant expansion is a mode 
sum, and may be carried out in those cases where the wave equation for the field 
mediating the self-force admits a mode expansion of the solution. In particular, this 
method can be used to calculate the gravitational self-force for a particle of mass n 
orbiting a black hole of mass M to order /i^, provided /i/M<Cl. We discuss how to 
use these two expansions to construct a full self-force, and in particular investigate 
criteria for matching the two expansions. As with all methods of computing self-forces 
for particles moving in black hole spacetimes, one encounters considerable technical 
difficulty in applying this method; nevertheless, it appears that the convergence of 
each series is good enough that a practical implementation may be plausible. 
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1. Introduction 

There is currently considerable interest in calculating the self-force on a charged particle 
moving in a curved background spacetime. This interest is twofold. On one hand, 
self-force is an intrinsically interesting phenomenon, apparently both deep and subtle. 
On the other, there is increasing practical interest in understanding the motion of 
particles with small but non-negligible masses in curved geometries. In particular, there 
is wide-spread belief that accurate computation of gravitational self-force is needed 
to calculate gravitational waveform templates for key data analysis efforts associated 
with interferometric gravitational wave detectors. With such templates, the Laser 
Interferometer Space Antenna (LISA) fll is expected to produce a wealth of information 
about strongly- field gravity jSj. 

We define self-force to be any force on a particle of quadratic (or higher) order in the 
charge carried by that particle. For electrically charged particles in fiat spacetime, the 
self- force is given by the Abraham-Lorentz-Dirac formula [SI EI- DeWitt and BrehmelH] 
derived the general expression for the self-force on an electrically charged particle in a 
curved background (however, they missed a term which was later supplied by HobbsjHI). 
The gravitational self-force was first calculated almost simultaneously by Mino, Sasaki 
and Tanaka [7j and by Quinn and Wald [8 . Later, Quinn derived the equivalent 
formula for a charge coupled to a minimally-coupled massless scalar field These 
results have resolved many of the issues of principle in computing self-forces in curved 
spacetime. PoissonfTU] has recently written a comprehensive review article in which 
these expressions and their implications are discussed at length. 

The self-force expressions for all fields of physical interest are of similar form. We 
henceforth restrict our attention, therefore, to the simplest particle with mass 

fi and scalar charge q coupled to a minimally-coupled massless scalar field in a curved 
background geometry. Lessons learned here should be extendible to other physical fields 
without major modification. 

The field equation for a minimally-coupled massless scalar field is 



Here, □ is the D'Alembertian of the curved background and p is the charge density. We 
consider a point particle, in which case 



where 5^(-) is a generalized Dirac distribution in four dimensions and z{t) denotes the 
worldline of the particle. 

For such a particle, Quinn^9j has shown that the self-force is given by 



□0 = -Anp. 



(1) 




(2) 




a 
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(3) 
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The quantities in equation are defined as follows: is the four-velocity of the 
particle and a° = u^V/jm" its four- acceleration. d° = w^V^a" denotes derivative of 
the acceleration with respect to its proper time, and = a°'aa is the magnitude of the 
acceleration squared. The quantity Rafs is the Ricci tensor of the background spacetime, 
and R is its scalar curvature, r is the proper time of the particle at its current position, 
while t' denotes the proper time at any other point along the particles worldline. Finally, 
Gj-et{x, x') is the retarded Green's function for the scalar field equation (fTJ), which satisfies 
the Green's function equation 



and has support only when x' is in the causal past of x. 

Notice that the terms on the right-hand-side of (jHl) can be gathered into three 
groups. The first group are local terms involving the acceleration of the particle. These 
give the Abraham-Lorentz-Dirac force on an accelerating particle in flat spacetime [Sj- 
They arise because the accelerating particle radiates, and that radiation produces a 
"radiation reaction" or recoil force on the particle. It is interesting to note that for a 
particle in Minkowski spacetime with constant a (whose worldline is a hyperbola), that 
d" 7^ (because the direction of the four- vector changes). In fact, a" = a^u" in this 
case, and the first group vanishes. This fact, which might seem surprising at first, can 
be understood as a consequence of the equivalence principle [TTl|T2llin]- We further note 
that for a particle freely falling in curved space (i.e. following a geodesic), a° = 0, and 
thus these terms will also vanish. 

The second group are also local terms, this time involving the background curvature. 
These terms involve only the Ricci curvature, and thus represent a self-force mediated 
by the matter content of the background. Interestingly, it is not necessary for the matter 
to interact directly with the particle for this to be true. We will largely ignore these first 
two groups of terms for the remainder of our discussion because for a particle following a 
geodesic in a vacuum background spacetime, which is the case of most practical interest, 
they vanish. Moreover, even when they do not vanish, they are easily calculated. 

In contrast, there is considerable practical difficulty in calculating the single non- 
local term which constitutes the third group. Interestingly, this non-local self-interaction 
arises in curved backgrounds, but not in a flat background. This is because of two ways 
in which the propagation of massless fields differs in curved and fiat backgrounds. In 
Minkowski space, massless fields propagate along null geodesies. Thus, a particle would 
have to be null-separated from some point on its past worldline to affect itself. However, 
the simple causal structure of Minkowski space does not allow this for a massive particle. 
The particle is restricted to a time-like geodesic, and no two points in Minkowski space 
can be connected both by a time-like geodesic and a null geodesic. Thus, any point with 
which the charged particle can interact, it cannot travel to, and vice versa. This is not 
true in a curved background, however, where the causal structure can be considerably 
more complicated. In this case, it is possible for the field, which "leaves" the particle 
along a null geodesic, to re-intersect that particle, which follows a timelike geodesic, at 




(4) 
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a later time. 

If this were the only mechanism by which non-local self-interactions arose, then 
there would be no self-interactions arising from within the normal neighbourhood of the 
point at which the particle sits. Recall that the normal neighbourhood of a point x is 
the set of all other points which are connected to x via a unique geodesic. Recall also 
that such a neighbourhood is guaranteed to exist. Thus, if a particle at x is massive 
and non-accelerating (i.e. following a timelike geodesic), then any point on the particle's 
worldline connected to x by a null geodesic is connected to x by (at least) two geodesies. 
Therefore, by definition, it is not in the normal neighbourhood of x. 

However, in general, fields do not propagate only along null geodesies in curved 
backgrounds. This is apparent from the Hadamard[T^ form of the retarded Green's 
function, 

Gret{x,x') = e[t-t']|f/(x,x') S[ct{x,x')]-V{x,x') 6 [-(t(x, x')] } , (5) 

where 9[-] is the Heaviside step function, 6[-] is the Dirac delta distribution, a{x,x') is 
one half of the square of the geodesic distance between points x and x', and U {x, x') 
and V{x, x') are smooth functions which depend on the details of the background and 
field equation. Note that this expression is only well defined if one can unambiguously 
define geodesic distance between x and x', which implies that x' is within the normal 
neighbourhood of x. 

Recall that the Green's function gives the field at position x due to a charge at 
position x'. Now, the first term on the right-hand-side of equation (jS)), which is known 
as the direct part of the Green's function, has support only when the geodesic distance 
between x and x' vanishes, or, in other words, when x and x' are separated by a null 
geodesic. This term is always present, but, as described above, cannot contribute to 
the self-force. However, the second term, which is known as the tail part of the Green's 
function, does contribute whenever the square of the geodesic distance between x and 
x' is negative, or, in other words, when x and x' are separated by a timelike geodesic:}:. 
Since every point on the particles past worldline is timelike separated from the particle 
by a timelike geodesic (to wit, the worldline), through this term a self-force can be 
generated at every time within the normal neighbourhood. 

Given that all points on a particles past worldline can interact with the particle, it 
may seem somewhat arbitrary to have singled out in the above discussion points on the 
past worldline which are null-separated from the particle. Indeed, in the literature, it is 
often stated that the self-force arises from the tail part of the field (or Green's function) 
and left at that. This statement can be somewhat confusing, however. Outside of the 
normal neighbourhood, there is no clear distinction between tail and direct parts - is 

I This fascinating fact was first elaborated upon by Hadamard^l|, who described it as a failure of 
Huygens' principle to hold in general for hyperbolic partial differential equations - c.f. the Paul Giinther 
memorial edition of Zeitschrift fiir Analysen und ihre Anwendungen [TS) for a collection of articles 
addressing the proof of the modified Hadamard's conjecture, which states that the only spacetimes for 
which standard wave equations obey Huygen's principle are those conformal to Minkowski spacetime 
and one plane-wave family of spacetimes. 
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the field from a point connected by both timelike and null geodesies a direct field, a tail 
field, or both? 

Furthermore, there is reason to believe that these points play a qualitatively 
different role in the self-interaction. In particular, when doing the integral over the 
Green's function, distributions (i.e. (5-functions and step functions) will be encountered 
when the source point and field point are null separated. In effect, the particle can 
"feel" its own direct field "sent" from points in the past. This is depicted in figure ^ In 
the specific case of the 0[M] Green's function, such distributions do appear, and their 
contribution make up the entire self- force (see equation (j211))- 

Let us turn now to the matter of calculating the non-local part of the self-force, 

r = g2 lim r V^Cct (^(r), z(r')) dr'. (6) 

The key feature of this formula is that the integral is well behaved over the entire 
region of integration, i.e. it is integrable throughout its domain. The termination of 
the integral at r — e (as opposed to integrating all the way to r) is a required for this 
to be true, however, because the Green's function has an unintegrable singularity at 
the coincidence limit (e = 0). Nonetheless, this contribution is distributional, and the 
integrand is perfectly well behaved as long as we approach r' = r from the past. 

At first glance, evaluating equation (jH} seems like a relatively straightforward task 
- an approximate retarded Green's function can be calculated, for instance, using mode- 
sum techniques ^H]- However, in practice, we can only sum a finite number of modes and 
thereby obtain an approximate Green's function. Furthermore, as we will demonstrate 
in section El the number of modes needed to obtain a given accuracy grows without 
bound as epsilon approaches zero. This does not preclude the use of modes to calculate 
the non-local part of the self-force, and, as can be seen in the pages of this special issue, 
modes are indeed widely used. The modes must, however, first be regularized by some 
method. 

Nonetheless, it is our purpose in this paper to explore an alternative method which 
may have some advantages over a regularized mode sum. It does not need to make use 
of regularized modes, although one might choose to do so. Rather, it is a method of 
matched expansions, in which one calculates the self-force using the tail within some 
portion of the normal neighbourhood of the particle and using an unregularized mode 
expansion for the remainder of the particle's worldline. More precisely, we propose to 
express equation ^ as 

r = -q' f V" V (^(r), z{r')) dr' + f ^V" G^et (^(r), z{r')) dr, (7) 

where Ar is an interval of proper time. We require that Ar be chosen so that z{t') 
is within the normal neighbourhood of z{t) for all r — Ar < r' < r. In other words, 
Ar distinguishes the contribution to the self-force coming from the recent history of the 
particle from the contribution that comes from the more distant past. It is the choice 
of Ar and the feasibility of evaluating the two integrals that will occupy us for the 
remainder of this paper. 
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Figure 1. In order to compute the force at z{t), one needs to integrate bacii over tiie 
entire liistory of the particle. Figure ^a depicts the basic idea behind this paper: use 
the Hadamard expansion to integrate over the portion of the worldhne "near" z(r), 
and then finish off the integral using a mode sum for points well separated from (and 
extending beyond the normal neighbourhood of) z{t). The choice of when to do the 
switch is also alluded to in the figure: the boundary of the normal neighbourhood - 
the crossing of geodesies - will likely be due to null geodesies looping around the black 
hole and re-intersecting the particle's worldline. This is discussed in Section [3 and 
the results are given in Table ^ It is probably good to do the switch somewhere in 
the middle, where keeping only a few terms in the Hadamard expansion is still a good 
approximation and before the mode sum is fighting against the divergent quantities 
on the light cone. FigureHb shows that for the 0[M] Green's function (see sectionEJ 
it is actually the portion of the past trajectory that is well removed from the current 
position (i.e. t' < t — [r + r')) that actually produces the "tail" force. Rigorously, 
it make little sense to talk of a "normal neighbourhood" for a spacetime with weak 
curvature, yet the similarity in figure ^a and figure strongly suggests that the 
portion of the worldline near the normal neighbourhood boundary and outside the 
normal neighbourhood boundary dominate the contribution to the self force. 

There are several notable features of equation (0). First, this expression is exactly 
equivalent to equation (jH)), since, as is evident from equation (jSJ, Gj-ct{x, x')= — V{x, x') 
provided we restrict ourselves to the interior of the past light cone, as is always the case 
for this integral. Second, there is no longer a limit needed because V{x,x') is regular 
everywhere. Finally, notice that we now have a new parameter, Ar, which we are free 
to choose so long as it is not too large. 

One might suspect that the contribution to the Green's function "falls off' fast 
enough, in general, so that perhaps only the first integral needs to be evaluated; and, 
since it is restricted to the normal neighbourhood, only the Hadamard expansion is 
needed to compute the self-force. Surprisingly, there are simple situations in which this 
line of reasoning turns out to be false; the second integral gives a significant (perhaps. 
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the dominant) contribution to the total force. We demonstrate this in section |21 using 
a Green's function for spacetime with a large central mass M where we only keep the 
terms of leading order in the central mass (the "0[M] Green's function"). Using the 
0[M] Green's function, none of the force originates from the tail contribution in the 
normal neighbourhood; all of the force comes times prior to the light reflection time. 

Unfortunately, neither integrand in equation ^ can be calculated exactly. As 
mentioned above, Gy-et{x,x') can be approximated in the second integrand by a mode 
sum expansion. The advantage here is that the mode sum does not need to be extended 
to the limit of the particle's position. Thus, for any fixed finite precision required, a 
finite number of modes will be needed for the second integral in equation (0). As noted 
above, however, that number will grow as At decreases. 

On the other hand, for V{x, x'), we can take advantage of the fact that we are within 
the normal neighbourhood, and can therefore define a Riemann normal coordinate 
system. In such a coordinate system, it is relatively straightforward to expand V{x,x') 
in a covariant Taylor series. Such a series will have coefficients constructed of geometric 
quantities (notably, the particle's four-velocity and the curvature of the background) 
evaluated at the particles position, and will be an expansion in geodesic distance from 
that position. Again, only a finite number of terms will be available, so only an 
approximation of V{x,x') can be calculated. In contrast to the mode sum expansion of 
G,:ct{x,x'), the covariant Taylor series will, in general, require more terms to achieve a 
given accuracy as At increases. 

The goal, then, is to calculate both series to sufficient accuracy within their 
applicable domains, and to then integrate and sum. We may take advantage or our 
parameter. At, to adjust the number of terms required in each series to achieve the 
required accuracy. However, both series are non-trivial to calculate, and the difficulty 
of calculating each successive term is greater than for the last. There is, therefore, no 
guarantee that there is any value of At such that the number of terms needed in both 
series can be calculated in practice. 

It is the main goal of the remainder of this paper to investigate if it is plausible 
that there would exist a choice of At which would make this scheme, first proposed by 
Poisson and Wiseman jl7j. viable for practical calculations. Because of the difficulty of 
performing both expansions, we will exploit existing results when possible. Furthermore, 
we will use the simplest and most transparent results that still bear upon the problem. 
We will not explicitly try to calculate the self-force for any specific geometry or particle 
motion, but will, rather, simply investigate the convergence of both series expansions. 

The plan of the remainder of this paper is as follows: we will review and examine 
the existing results for the expansion of the V{x, x') in the next section. Following 
that, in sectional we will investigate the mode expansion. Finally, we will discuss our 
conclusions in section |3] 
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2. The quasi-local part of the self- force 

Covariant normal neighbourhood expansions have a venerable history for calculations 
of both self-forces IS] and quantum field effects in curved backgrounds [THl IT^ . There is 
therefore a plethora of material from which to begin a calculation of the first term in 
equation (jZj), which we will call the quasi- local part of the self- force 



/ql = -?'/ V-V{z{r),z{T'))dT'. (8) 



'ql 

J r-Ar 

There are, however, to our knowledge, only three papers in the literature which address 
the calculation of this part of the self-force. The first was by Roberts, in which he 
stopped just short of calculating the quasi-local part of the self-force for an electric 
charge in an arbitrary background |2nj. A trivial extension of his paper gives the leading 
order to the quasi-local part of the self-force in this case, 

/ql = - + ucu^) C%,',su^u'< At' + 0{At'). (9) 

Here, Ca/s^s is the Weyl curvature tensor of the background at the particle's position, 
is the particle's four-velocity, and At, as previously, is a proper time interval along 
the particle's past worldline. 

The second result is due Anderson and Huj21j. They calculated the tail of the 
retarded Green's function in the normal neighbourhood for a particle with a minimally- 
coupled massless scalar field in a Schwarzschild background geometry. By using a 
Hadamard-WKB expansion for the Euclidean Green's function, they are able to expand 
V{x, x') to sixth order in the geodesic separation. They do not go on to calculate the 
quasi- local self-force explicitly (although that is clearly the motivation for their paper), 
but it is a relatively straightforward matter to do so. Because their results are somewhat 
unwieldy, we refer the reader to their paper for the actual result, although we will be 
reproducing the self-force derived from it for two simple particle motions below. 

Most recently, the quasi-local part of the gravitational self-force for a particle 
in an arbitrary curved background has been calculated by Anderson, Flanagan and 
Ottewill|22j. They calculate the first two non- vanishing terms in the Taylor series, and 
find 

fQLair,AT) = -ii\5l^ + u^u^)Cp^seC,\'u'u''uPAT' 



1 3 



19 /I 



Ar-^ + 0(Ar*). (10) 



60 V2 

They also provide explicit expressions for some particle motions in Schwarzschild and 
Kerr geometries. 

For the purposes of evaluating the convergence of these expansions, we will focus 
on the results of Anderson and Hu||21|, both because the scalar field case is inherently 
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the simplest and because the provide the highest order expansion with which to work. 
Unfortunately, Anderson and Hu do not take the derivative of V{x, x') necessary to 
calculate the self-force contribution. Furthermore, their result is expressed in terms 
of coordinate expansion rather than an expansion in geodesic distance. It is only for 
relatively simple particle motions that one can easily recast coordinate distance into 
geodesic distance. Fortunately, it is in the spirit of our explorations here to take such 
simple cases, and we shall. 

The fundamental question we ask in this section, then, is "what is the rate of 
convergence of the expansion for the quasi-local part the self-force as a function of At?" 
Clearly, as Ar— >0, the expansion becomes exact at any order. The more interesting limit 
is the one where At approaches the boundary of the normal neighbourhood. It might 
be reasonable to expect that the expansion ceases to converge at all in that case, since 
that is the boundary of its domain of validity. Indeed, there is no fundamental reason 
that the domain of convergence of the series expansion could not be much smaller than 
the normal neighbourhood. Addressing this question, however, is slightly complicated 
by the fact that it is not immediately clear what the value of At at the boundary of 
the normal neighbourhood is. 

In the case of a Schwarzschild black hole spacetime, at least, one can get some 
insight into the normal neighbourhood boundary by studying null geodesies intersecting 
a particle in a circular geodesic at radius R. For such a particle, the angular displacement 
is related to the particle's proper time by 



_ , M At 

A0particle - y -^^^ ^, (11) 

where M is the mass of the black hole. On the other hand, the angular deflection (as 
seen at r = i?) of a null geodesic passing within a nearest distance Rq > 3M of the 
black hole is 



A0p.oton = 2 J du ^ ^_^M-Riu\l-2Mur ^^'^ 

(note that if Rq < 3M, there is no turning point for the geodesic and it does not return 
to larger radii). Finally, the proper time, as measured by the particle, for a photon to 
descend from a distance R to within a distance Rq of the black hole and return to a 
distance R is given by 



Ro-2M\ r I r^{Ro-2M) 



H ^0 ) J no'^\-2M\JrRoir^-R'o)-2M{r^-Riy ^^^^ 

What we would like to find is the minimum value of Rq (and hence At) for which a null 
geodesic intersecting the particle's circular geodesic at a given R can re-intersect it at a 
the same R but a different time. For a particle orbiting at i?^3M from the black hole, 
the condition for this to occur is 

A^paxticle + A^photon = StT. (14) 
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R/M 



Ro/M 



At/M 



6 

10 

20 
100 
1000 



3.46 
3.55 
3.58 
3.56 
3.54 



18.6 
31.9 
57.8 
228 
2038 



Table 1. Values for a null geodesic leaving a particle in a circular geodesic orbit in 
Schwarzschild spacetime at radius R and returning to the particle. Rq is the minimal 
radial coordinate value of the null geodesic (minimum coordinated distance from the 
centre of the black hole) and Ar is the proper time the particle has to wait for the null 
geodesic to return. 



Assuming this and substituting equations (fTT] - fT!?|l into equation (|T^ . we can choose a 
value of R and solve equation (fT^ numerically for _Ro- We give some values obtained in 
this way in the table [H 

There are a number of noteworthy features of tabled First, at at approximately 
Ro = 3.2M, the angular deflection of the null geodesic is 27r, and this therefore represents 
the value of Rq corresponding to a particle orbit at R-^oo. Next, we note that the values 
of Ar take values between ~3i? for the closest orbits and ~2i? for the most distance 
orbits, asymptoting to Ar = 2R in the limit R—>-oo. In other words, for orbits at large 
distances the time is dominated by the time for the photon to reach the black hole and 
return. 

The Ar's quoted in table represent approximate upper bounds on the extent of 
the normal neighbourhood along the past worldline of the particle since they demarcate 
two intersections of the particle's geodesic with a null geodesic. Thus, for a particle 
at a fixed radius R, we need not worry about the convergence of the expansion of the 
quasi- local part of the self- force beyond ~2i? to the past. 

Let us now consider two such expansions. First, we consider the expansion for a 
static particle at radius R coupled to a minimally-coupled massless scalar field by scalar 
charge q. We take the expression for V{x,x') given by Anderson and Hu|2I] and take 
partial derivatives with respect to the Schwarzschild coordinates. Next, we convert these 
coordinates into proper time using the coordinate parameterization 



which is appropriate for a static particle. We also, without loss of generality, set ^ = 7r/2. 
We can then integrate the expansion of VaV{x, x') with respect to proper time r to 
obtain 



Act) = 0, A^ = 0, Ar = 




(15) 




.7 



(16) 
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which is the only non-vanishing component of the quasi-local part of the self-force to 
the order of this expansion. 

Denote the fifth and seventh order terms in equation (jlfjj) as 

/ql[5]=^^ m - IIM) {R - 2Mf Ar^ (17) 
1 fl^M^ 

/ql[7] = oo^^ %..n (20i?^ - 195Mi?2 + 598M^R - 585M^) {R - 2Mf At\ (18) 
3360 R'^^ 

Then the fractional truncation error which is induced by only taking terms in the series 
expansion of the quasi-local part of the self-force to order Ar^ can be estimated as 

/ql[7] 



^"/ql[5] + /^l[7]- ^^^^ 

This gives an estimate of the upper limit on the local truncation error (the error in 
truncating the next term in the series) rather than the more desirable global truncation 
error (the error in truncating all remaining terms in the series). Nonetheless, it is a 
standard measure of truncation error for this kind of analysis where neither the exact 
solution nor a form for the general term in the series is known, and is in any case the 
best estimate of truncation error we have for this series. 

As noted previously, the truncation error can be expected to grow with Ar. We 
have calculated the error as a function of Ar for particles located at r = 6M, lOM, 
20M and lOOM. The results are presented in figure El We see that, despite our concern 
that the quasi-local expansion might fail to converge at the boundary of the normal 
neighbourhood, it does, in fact, appear to converge everywhere within the normal 
neighbourhood. Furthermore, it converges quite well for particles at all radii, with 
estimated fractional error less than 0.1, back to almost Ar = lOM. For R = QM, which 
is the closest particle that we consider and also the case for which the self-force should 
be most important in calculating templates for LISA, the fractional error is less than 
0.1 more than half way out to the boundary of the normal neighbourhood. 

We can apply exactly the same method to a particle in a circular geodesic orbit at 
radius R around Schwarzschild. In this case, the coordinates are related to proper time 
by 



Ae = 0, Ar = 0, A0 = -W — ^^(r-r'), At = J (r-r'). (20) 

^ R\ R-3M^ \ R-3M^ ' ^ ' 

Again, we have set = 7r/2. In this case, there are three non- vanishing components of 
the quasi-local part of the self-force 



, _ g^M^ / R 1 
~ i?i2 \l R-3M{R-3My 



X 



^ R^{R-2M){R-3M){5R-19 M) Ar^ 



2240 



(56R^ - 728 R^M + 3461 R^M^ - 6990 M^R + 5073 M^) Ar^ 
13440 ^ ^ 
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Figure 2. Estimated fractional error e as a function of At in the radial component 
of the self-force on a static charge coupled to a massless minimally-coupled scalar field 
in a Schwarzschild spacetime. We show a particle at four different distances from the 
black hole. Upper bounds for the range of At are those calculated in tabled 
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Estimated Error in Quasi-Local Self-Fojce 

Time Component for Circular Orbit at Order At 




At (M) 



Figure 3. Estimated fractional error as a function of At in the time component of the 
self-force on a charge coupled to a massless minimally-coupled scalar field orbiting on 
a circular geodesic in a Schwarzschild spacetime. We show a particle at four different 
distances from the black hole. Upper bounds for the range at each distance are again 
approximately those from table ^ 



(23) 



We compute the estimated fractional error for each of these components in the 
same manner as we did for the static particle. The results are presented in figures El IH 
andEl Again, we see convergence everywhere within the normal neighbourhood for all 
components of the quasi-local part of the self-force, and good convergence up to fairly 
high values of At. 



3. The contribution to the self-force from the distant past 

In the previous section, we examined the contribution to the self-force from the portion 
of the worldline within the normal neighbourhood of the field point at z{t), i.e. from the 
first integral in equation ((7j). In this section, we will examine the contribution from the 
earlier part of the trajectory. In particular, we will examine the feasibility of computing 
the second integral in equation using a mode sum expansion for the Green's function. 
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Figure 4. Estimated fractional error as a function of At in the radial component of 
the self-force on a charge coupled to a massless minimally-coupled scalar field orbiting 
on a circular geodesic in a Schwarzschild spacetime. We show a particle at four different 
distances from the black hole. 



3.1. The 0[M] Green's function 

In studying the forces on a freely falling electric charge in a curved spacetime, DeWitt 
and DeWitt [22] developed some very clever techniques for finding the Green's function 
for a spherically symmetric spacetime with mass M at the centre. They give an 
approximate expression - accurate to leading order in the central mass - for the Green's 
function for an electric charge in this Schwarzschild-like spacetime. This method was 
later extended by Wiseman [21] and Pfenning and Poisson [23] to the case of a scalar 
charge. The Green's function is found to be 



X — X' 



+M— L ^[t-t'-\^-^'\] In r r + r' + {t-t') \ _^m- t') -ir + r' 



'12 



dt' { |x-x'| \r + r' - {t -t') J (t - t')^ - |x - x' 

+0[M^] , (24) 

where (t,x) is the field-point and (t',x') is the source point. The first term is clearly 
identifiable as the retarded Green's function in Minkowski (fiat) spacetime. The 
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Estimated Error in Quasi-Local Self-Fprce 

(|) Component for Circular Orbit at Order At 




At (M) 



Figure 5. Estimated fractional error as a function of At in the (p component of the 
self-force on a charge coupled to a massless minimally-coupled scalar field orbiting on 
a circular geodesic in a Schwarzschild spacetime. We show a particle at four different 
distances from the black hole. ^ 



logarithmic term is a correction to the light cones§. The final term is the "tail" of 
the Green's function. Notice that this term only contributes prior to t' = t — {r + r'), 
the light reflection time. This is depicted in figure H This gives a strong indication that 
the self-force is dependent on the portion of the worldline that is outside the normal 
neighbourhood. 

Multiplying by spherical harmonics and integrating over the solid angle, we can 
obtain the angular mode decomposition of this Green's function 

xP;(0 e[t - t' - |r - r I ] 0[r + r' - {t - t') ] 
+AM^Qi{0 e[t + t' - (r + r')]| P^[cos(7)] 

+ 0[M^] , (25) 

§ This term reflects the fact that the light cones are bent in curved spacetime. It was not included in 
the DeWitt-DeWitt EBl calculation. 



dt' \r + r'-{t- V) J 
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Figure 6. The point cliarge only lias an instant of causal contact with the field 
point {t, r) when the past light cone intersects the worldline of the charge; however 
when using multipoles, the point charge is smeared over a sphere and therefore has an 
extended period of interaction. In equation H29|l we have placed the field point on the 
worldline of the source, i.e. 7 = and r = b. In figureElwe show how the gradient of the 
field "builds up" for a partial sum of multipoles as we integrate over t' . In computing 
the numbers in table [3 we have again placed the field point on the worldline of the 
source and we have integrated from (3 — n to (3 = tt/2, i.e. cos/3e[— 1, 0]. 



where Pi and Qi represent Legendre functions and 

J, ' =-^°^g (26) 

COS 7 = cos 9 cosO' + sin sin 6' cos(0 — 0') (27) 

and (t, r, 9, (p) and {f, r', 6', cf)') are the spherical coordinates of the field-point and source- 
point respectively. The angle (3 is shown in figure El In equation (j2Sl), it is understood 
that the first partial derivative not only operates on the natural log, but also on Pi{^) 
and the 9-functions. 

In numerical implementations, we would only be able to include a finite number of 
multipoles. In this case, we can avoid doing a numerical summation by computing the 
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partial sum analytically, i.e. 



N 

G{N; t, r, t' , r', cos 7) = terms in equation 

1=0 




d_ r i + (iV + l)(P7v+i(cos7)g^(0-P7v+i(cos7)Q7v+i(0) Q.^ . ^.1] 

9t' [ ^-cos7 ^ ^ ) \ 

+0[M^] . (28) 



As above, the first partial derivative acts on the natural log as well as the Legendre 
functions and the B-functions. 

This formula can now be used in the second integral in equation 0. As in Section|21 
we focus our attention on a static source at radius R and cos 7 = 1 (see figure IHl). Notice 
that the final term - the tail term - is a total derivative. When this is substituted into 
equation ((7j), we can integrate by parts and see immediately that the tail contribution 
to the force vanishes. This is just a re- confirmation of the well known result that there 
is no self- force on a static scalar charge in Schwarzschild spacetime (c.f. Wiseman |26j). 
By the same argument, the natural log term also gives no contribution to the force 
integral. 

Substituting the remaining term - the flat-space portion of equation ()28|) - into 
equation (jZj), a straightforward calculation reveals that the first multipoles of the 
source give a radial component of the force 



where = 1 — M/b + 0[M^] is the leading order contribution to the time- component 
of the four-velocity of the static charge. The upper limit is defined as 



As can be seen from figure when the field point is on the worldline of the charge, cos (5 
can be used to parameterise the portion of the worldline from the field point (cos /3 = 1) 
backwards to the time when the past null cone of the field point first intersects the 
cylinder around the central mass on which the worldline lies (cos/3 = —1). As we will 
discuss below, this is the only portion of the worldline which can contribute to the self- 
force in this case. In our case, we do not integrate to the field point, but rather to some 
time At to the past of it as per figure ^ thus terminating our the integral at cos/5max- 
This is indicated on the left hand side of equation (j^^ by the argument At. 

The only portion of the Green's function that contributed to equation (jSHl) was the 
fiat-space term. Since the static point particle would never intersect the past light cone 
of a field point located on its own worldline, we can ask: why doesn't /^(Ar, A^) vanish 




(29) 




(30) 
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identically? The answer lies in the fact that we are using only a finite number of modes 
to describe the field gradient. 

The decomposition of the field into angular modes can also be thought of as a mode 
decomposition of the source. In the case of the static point particle, we are, in effect, 
replacing the point charge with a sum of spherical shells of charge where each shell has 
radius R and an angular distribution of charge smeared on it. The point particle is only 
recovered when a large (infinite) number of shells of charge are summed. As is evident 
from figure El although the point charge doesn't come into null contact with the source 
point, the spherical shells of charge do intersect the past light cone of the field point. 
If we terminate the integral at some At > 0, but include an infinite number of terms 
in the summation (i.e. N—>-oo), we will recover the point-particle nature of the source, 
and therefore there would be no contribution to the force. However, in any practical 
calculation, we will need to terminate at some finite value of A^, and thus we will be left 
with an unwanted, and inescapable, contribution to the force. The key question is can 
we include enough terms in the summation to squeeze this unwanted contribution to an 
acceptably low level? 

Figure U\ shows how the field gradient accumulates as we integrate over r' for a 
fixed number of modes. The oscillations are artefacts of the finite number of modes 
summed, and it is the envelope of the function which reflects the level of accuracy that 
can be achieved. If we are to accurately evaluate the value for the self-force, we must 
squeeze the envelope of this function below the value of the self-force. For example, on 
dimensional grounds, one might expect there to be a radial component of the self-force 
for a static scalar charge in Schwarzschild spacetime of the form 

= ^ V ' ^^^^ 
where A is an unknown coefficient we would be trying to dig out from under the unwanted 
contribution to the force from our finite mode sum (in the electrostatic case, there is 
a finite contribution to the self-force of exactly the form shown in equation (pTT|) with 
A = 1). The amplitude of the envelope is therefore an indication of the bound we place 
on the contribution to the self-force. In the present case, we know the solution vanishes 
identically, and evaluating equation (fTT)|l from section |21 at i? = 6M, it can be seen 
that the first integral in equation ((7j) gives only a tiny contribution to the self-force of 
a static scalar charge, even when At is extended out to near the edge of the normal 
neighbourhood. Therefore, in this case, the amplitude of the envelope is nothing but 
the approximate accuracy to which we have calculated the self-force. 

If we choose At such that cos/5max = 0, which is about half-way out to our upper 
bound on the edge of the normal neighbourhood, then the we see from table 121 that 
using 200 multipoles we can constrain the A to be less than 0.1. 
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Figure 7. This sliows ttie accumulation of the field gradient as we integrate over 
t'. The normalization is {APu* /q'^)fr{AT, N). Here N is chosen to be 50. If we had 
chosen a larger N, the frequency of oscillation in the plot would increase, but the 
envelope of the oscillation would become smaller. The numbers in the table El are a 
measure of how fast the envelope shrinks as the number of modes is increased. For 
example the entry for 50 modes is just the value of the peak near cos /3max — 0. 



Modes 


(MV/g2)/,(Ar,iV) 


bound on A 


10 


0.0017 


0.36 


50 


0.00078 


0.17 


100 


0.00056 


0.12 


200 


0.00041 


0.088 


300 


0.00036 


0.078 



Table 2. Table showing the convergence of the mode sum and the accuracy to which 
we could use it to compute the gradient appearing in the integrand of equation Q. The 
first column is the number of modes included in the sum. The integral is terminated 
at At chosen such that /?max = 0. The charge is located at i? = 6M. This value in 
the second column is the peak nearest cos /3max = on plots similar to figure [7| The 
right column shows the accuracy with which you could constrain the coefficient of a 
self-force of the form \q^M/b^ with the number of modes used. 
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4. Conclusions 

Our goal in this paper was to assess the feasibility of the Poisson- Wiseman matched 
expansion scheme for calculating the self-force. This scheme involves splitting the non- 
local term in the self- force (which is an integral over the particles worldline z{t)) at 
some proper time r — At. For the — oo < r' < r — At (distant) part of the integral, 
the integrand is expanded as a mode sum. For the r — At < t' < t (quasi-local) part, 
the integrand is expanded in a covariant Taylor series, which requires z{t — At) be in 
the normal neighbourhood of z{t). For the resulting self-force to be accurate, the split 
must be done in such a way that each expansion reach sufficient accuracy with a finite 
and feasible number of terms. Until now, it has not been clear that such a split even 
existed. 

For this preliminary investigation, we have studied the simplest scenarios: 
minimally-coupled massless scalar field, Schwarzschild geometry, static particles and 
circular geodesic motion. Our results are somewhat surprising. The quasi-local 
expansion appears to converge well - when the particle is at 6M and expanding to order 
Ar^, we achieve an estimated truncation error of the order of a few percent for Ar~6M. 
This would naively have seemed to us enough of a buffer between the mode-sum integral 
and the singularity to allow rapid convergence of the mode sum as well. 

However, we have found the convergence of the mode sum to be poor. This is 
because the mode sum smears the charge of the particle over spheres of finite radius 
which extend outside of the normal neighbourhood. Thus, even the fiat-space term 
in the Green's function, which cannot contribute to the self-force for a scalar particle, 
seems to contribute at every order. For our simple case of a static particle at 6M and 
Ar~6M, going from 10 modes to 100 increased our accuracy by only a factor of 3. 

Clearly, we need a way to speed convergence of mode sum. In the case presented, 
one could have achieved infinite accuracy at every order by regularizing the modes, 
which would have removed the fiat-space portion of the Green's function. It might, 
therefore, be possible to speed the convergence of the mode sum by regularizing the 
modesj^Z]. This would be ironic, since this method is supposed to provide an alternative 
to mode regularization. Nonetheless, it might be that combining regularization with 
a two expansion approach will provide much better convergence than either alone. 
Alternatively, we note in figure [7| that the mode-sum integrand oscillates about its 
true value. This leads us to speculate that averaging over these oscillations would lead 
to much quicker convergence of the mode-sum integral. 

While we believe our convergence analysis is general, and will apply beyond the 
simple scenarios we have explored, we would be remiss in not pointing out that there 
are additional issues for higher spin fields. The foremost of these is the issue of 
gauge - in order to meaningfully combine the quasi-local and more distant self-force 
contributions, they will need to be expressed in the same gauge. For instance, the 
quasi-local contribution for the gravitational field |22j has only be derived in the Lorentz 
gauge, where no expression for the modes is available. We note, however, that this is a 
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generic difficulty for gravitational self-force calculations because the formal expressions 
for the self-force are themselves derived in the Lorentz gauge [3 |H], and we are heartened 
by current progress in understanding gauge transformations in this context |2H]- 

In any case, we feel that these preliminary results are promising enough to warrant 
further investigation. A good figure of merit for calculational schemes like this is 
accuracy per floating point operation. As is, this approach would seem to lag approaches 
like mode-sum regularization when measured on this scale. Nonetheless, it might, for 
instance, provide con&mation for results from other approaches. Further, as mentioned 
above, it is still possible that this calculational scheme, with refinements, could rival or 
exceed in computational efficiency other known schemes, especially since the quasi-local 
expansion, once calculated, can be applied with little further effort to any spacetime 
or particle motion. In the mean time, this approach remains, in our opinion, in the 
category of "promising but possessing some technical challenges" . 
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